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ABSTRACT 

Planets are known to orbit giant stars, yet there is a shortage of planets 
orbiting within ~0.5 AU (P< 100 days). First-ascent giants have not expanded 
enough to engulf such planets, but tidal forces can bring planets to the surface of 
the star far beyond the stellar radius. So the question remains: are tidal forces 
strong enough in these stars to engulf all the missing planets? We describe a 
high-cadence observational program to obtain precise radial velocities of bright 
giants from Weihai Observatory of Shandong University. We present data on 
the planet host Beta Gem (HD 62509), conhrming our ability to derive accurate 
and precise velocities; our data achieve an rms of 7.3 ms“^ about the Keplerian 
orbit £t. This planet-search programme currently receives ~100 nights per year, 
allowing us to aggressively pursue short-period planets to determine whether they 
are truly absent. 

Subject headings: planetary systems - techniques: radial velocities - stars: giants 


1. Introduction 

Searches for planets orbiting evolved stars have been underway for more than a decade. 
A major science goal of these campaigns has been to explore the dependence of planetary sys- 
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tern properties on host star mass - precise radial velocities can be obtained for these “retired 
A stars” whereas A-type dwarfs present few absorption lines from which Doppler information 
can be extracted. However, of the 94 planets known to orbit giant stars (log g < 3.5), only 
seven have semimajor axes a < 0.5 AU. This is in spite of a strong bias favoring the detection 
of short-period planets. Abont 2000 evolved high-mass stars are currently being monitored 


worldwide by various teams, ( 

Hat 

zes & Goehr an 

1993 

Erin 

c et al. 

2002; 

Sato et al. 

2003; 

Setiawan et al. 

2003; 

Hatzes et al. 

2005; 

Johnson et al. 

2006; 

Niedzielski & Wolszczan 

2008; 

Lee et al. 

2011 

Wittenmver et al. 

2011bl. e.g.). Given that short-period planets are com- 


mon around solar-mass stars, there is an open question: Where are the short-period planets 
around high-mass stars? 


Villaver &: Liviol (120091) proposed that planetary orbits are affected by the evolution of 
the stars, showing how tidal interaction can lead to the engulfment of close-in planets. This 
process is strongly influenced by stellar evolution details such as the stellar mass, mass-loss, 
and metallicity. E ngulfment appears to be much more efficien t for more-massive planets and 
less-massive stars flKunitomo et al.ll201ll: IVillaver et al.ll2014f) . 


Giant and subgiant stars sample pl anet hosts that, in p r inciple, are more massive than 
their main sequence cou nterparts (e.g. 


Johnson et al. 

2010; 

Sato et al. 

2008; 

Bowler et al. 


2010 


Gettel et al.ll2012l) . This supposition (takei i alone) has been used to suggest a rela¬ 
tion between planet formation and stellar mass flGurriel l2009h . We note that the masses 


of these evolved star s have be en the subjec t of some controversy, e.g. iLlovdl (120111.120131) : 


Johnson et al.l (l2013l) . In brief, iLlovdl (1201 ll) argues that due to the stellar initial mass func¬ 


tion and the rapid evolution of more massive stars, a given evolve d star is more l ikely to be 
a “retired” Sun-like s tar than a “retired A star” as put forward by l.lohn.son et ahl (120061) and 


Johnson et al.l (120131 ). While this relation might still hold, there is, however, an increasing 


number of close-in planets being detected in the main-sequence stage orbiting A-F stars (e.g., 
HAT-P-49, Bieryla et ah 2014; WTS-lb, Gappetta et al. 2012; Kepler-14b, Buchhave et ah 
2011; WASP-33b, Gollier Gameron et al. 2010; KELT-3b, Pepper et al. 2013; OGLE2-TR- 
L9b, Snellen et al. 2009). These close-in planets are not found in radial velocity searches 
around evolved stars, suggesting that it is the evolution of the star, and not its mass, that 
plays a major role in removing planets from close orbits. Thus the discovery of a population 
of short period planets around evolved stars can offer fundamental insights into the strength 
of the tidal forces. Understanding these tidal forces is critical for producing accurate models 
of planetary orbits as stars evolve off the main sequence. Furthermore, they can help to 
shed some light on the dependency of the planet formation process on the stellar mass, with 
important implications for gas giant planet formation and migration. 


In fact, close-in planets are expected to enter the stellar envelope as the star leaves the 
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main sequence and evolves onto the red giant branch fjVillaver fc Livid l2009l: IVillaver et ah 


2014J) . As the star evolves, it removes planets from a region that extends far beyond the 
stellar radius to the entire region of tidal influence (a/i?* ~ 2 — 3 for a Jupiter-mass planet). 
During the subgiant phase, the star hrst clears out the very close-in planets present during 
the main-sequence evolution, and then proceeds to clear out a larger region as the stellar 
radius increases when it ascends the red-giant branch. However, the vast majority of the 
evolved stars targeted by radial-velocity surveys have not completed their ascent up the red- 

le evolved stars currently known to 


giant branch, and so are not at their maximum radius. T 
host plane ts have typical radii srn aller than about bR© fISato et al.ll2010l) . Stellar evolution 


models bv IVillaver &: Liviol (120091) show that for a 2Mq star, the radius exceeds 0.1 AU (i.e. 
a 10-day orbital period) for only the 40 million years immediately preceding the helium flash. 
For these reasons, we only expect planets orbiting closer than a ~0.15 AU to be lost to thei r 
expanding host star, and this only excludes planets with P < 15 days (I Johnson et al.l 120071) . 
Our survey will search for the planets with 15 < P < 100 days which are common around 
main-sequence stars, but missing around evolved stars. 


This paper is organised as follows: Section 2 describes the Weihai Echelle Spectrograph 
(WES), and the target selection criteria. Section 3 demonstrates the ability of the WES 
to deliver precise radial velocities by presenting new data and orbital hts for the known 
planet-hosting giant Beta Gem (HD 62509). Finally, we give our conclusions in Section 4. 


2. Observational Program 
2.1. The Weihai Observatory Echelle Spectrograph 

The Weihai Echelle Spectrograph (WES) is a bench-mounted, stabilised, hbre-fed spec¬ 
trograph attached to the 1-metre telescope located at the Weihai Observatory of Shandong 
University in Weihai, China. The spectrograph is thermally stabilised to ±0.10 K. Light is 
fed from the telescope into the spectrograph by ~10m of circular hbre with diameter 70/im 
via a 160/rm pinhole. It can achieve a maximum resolution of 57,000 with 2.2-pixel sam¬ 
pling. The WES is the hrst hbre-fed echelle spectrograph in China, and has the primary 
function of radial-velocity planet searc h. Further details about the spectr ogra ph and the 


Weiha i Observatory site can be found in lCao fc RenI (j2014l) . lHu et al.l (120141) . and I Guo et ah 

(I2OI4I) . 


The observing and data analysis procedures are typical of precise radial-velocity planet 
searches. Doppler velocity observations are performed at a resolving power of R ~45,000. 
Calibration of the spectrograph point-spread function is achieved using an iodine absorption 
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cell temperature-controlled at 65.0±0.1°C. The iodine cell imprints a dense forest of narrow 
absorption lines on the stellar spectrum, from 5000 to 6200 A, permitting the contemporane¬ 


ous c alibration of the spectrograph point-spread fun ction fIValenti et al.l Il995l: iButler et ah 


19961) . Velocities are obtained using the Austral code flEndl et ahll2000l) . w hich has been suc¬ 


cessfully used by several planet-search programs for more than 10 years (lEndl et al.l I2004J: 
Wittenmver et al.N2011bl: iRobertson et al.N2012aj) . The iodine region is broken into ~380 


100-pixel chunks, corresponding to about 4.3 Aper chunk. The chunks are weighted by their 
Doppler content (dehned as the sum of all pixel to pixel gradients), and each chunk pro¬ 
duces a radial velocity relative to the iodine-free template spectrum. The hnal velocity is 
computed as the mean value of the chunks after an iterative 3a clipping to reject outliers. 
The uncertainty is the standard error of the mean, i.e. the rms scatter divided by the square 
root of the number of accepted chunks. 


2.2. Target selection and observing strategy 


We chose a small sample of bright giants, starting with Northern hern isphere Hipparcos 
stars with V < 5.0 and luminosity classes III and IV (Ivan LeeuwenI 120071) . After removing 
all stars flagged for variability, double/multiple stars, and suspected composite spectra, we 
imposed a colour cut 0.5 < (B — V) < 1 . 2 to match the selection cri t eria of many evolved-star 


planet-search surveys ( e.g. ISato et al.l 120051: l.lohnson et al.l 120061: iNiedzielski fc Wolszczan 


20081 : I Jones et al.l 120111) . We then required My > 0.0 to exclude the most-evolved stars 


which w ould be subject to large-scale pulsations that confound radial-velocity searches for 
planets (IHekker et al.l 120081 ). To reduce the impact of pulsations on the detectability of 
planets, we also required that the Hivvarcos photome t ric sc atter be smaller than 0.005 
mag. From the scaling relations of iKjeldsen fc Bedding fjl9951) . this threshold corresponds 
to velocity amplitudes smaller than 5-10 ms“^, a level comparable to our photon noise and 
small compared to the expected amplitudes of giant planets in short-period orbits. The 42 
remaining targets are enumerated in Table [H 

Our knowledge of short-period planets orbiting these evolved stars remains limited by 
the observing strategies employed by the major programs. Radial-velocity surveys are usually 
subject to the exigencies of telescope scheduling, such that they are typically allocated time 
in a single short run each month (during the bright lunation). The “Rocky Planet Search” 
campaigns conducted by the Anglo-Australian Planet Search addressed this problem by 
observing a subset of ~30 stars every night for 48 consecutive nights. Th ose campaigns 


have demonstrated dramatically increased sensit ivity to short-period planets flO ’Toole et ah 
20091: IVogt et al.l l2010l: IWittenmver et al.l l2011aj) . Clearly, high cadence is the way forward 
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flSwift et al.ll2015l) . For example, IJones et al.l (120151) were able to conclusively detect the 


89-day planet orbiting the K giant HD 121056 by densely samplin g its orbit using CHIR ON, 
a spectrograph dedicated to precise radial-velocity planet search flTokovinin et al.ll2013l) . To 
achieve the best possible cadence, our survey has been allocated ~100 nights per year on 
the Weihai Observatory Im telescope. 


3. Preliminary Results: Beta Gem’s Planetary Companion 


Six of our 42 targets are known planet hosts (Table [2]). They are useful 


the sample because (1) furthe r monitoring can often reveal additional planets fjWright et ah 


members of 


20091: IWittenmver et ahll2013[) . and (2) the repository of data on these systems can serve as 
a check on the new data being obtained by this survey. 

Beta Gem (HD 62509, HIP 37826) is an extremely bright KO giant {V = 1.14) known 


to host a planet with P = 590 days and m sin i=2.76 Mjup (IHatzes fc Cochran! Il993 


Hatzes et al.l 120061) . We have 38 observations of HD 62509 over a 500-day baseline, with 
a mean internal velocity uncertainty of 8.6ms“^. Exposure times ranged from 300-900 sec¬ 
onds, with typical S/N~200-300 per pixel. An iodine-free template spectrum was obtained 
on 2014 March 7, and all velocities given in Tableware computed relative to that template. 

To check the consistency of the Weihai data, we include all published velocities span¬ 
ning more than 25 years to fit a Keplerian orbit to the plan etary signal. Our an alysis uses 
eight data sets, seven of which have already been published. iLarson et al.l (119931) presented 
velocities from the CFHT {N = 39) and DAO (A^ = 27), calibrated using an HF gas cell 
(jCampbell fc Walkerl 119791) . a prec ursor to the currently-use d iodine cell technique. The 
Beta Gem planet discovery paper (IHatzes fc CochranI 119931) gave 38 radial-ve locity mea¬ 
surements from the McDonald Observatory 2.1m telescope. Hatzes et ah (j2006 ) reported a 
furt her 22 epochs from the ongoing McDonald Observatory 2.7m “Phase 3” planet-search 
(e.g. Robertson et ah 2012a 3) and 11 epochs using a higher-resolution, limited wavelength 
setting on the same telescope (“cs21” as desribed in Hatzes et ah 2006 and Wittenmyer et al. 


2006). The Tan 
as presented in 


enburg Obse r vatory Planet Search (TOPS) also yielded 22 radial velocities 


Hatzes et al.l (120061 ). Finally, 80 velocities from Lick Observatory given in 


Reffert et ahl (120061) were included in our fit. 


We used the GaussFit nonlinear least-squares code (Ijeffervs et al.l Il987|) to obtain a 
Keplerian model ht for all eight data sets simultaneously. The velocity data now span 
33.5 years. Uncertain ties were estimated using the bootstrap routine within Systemic 2 
(iMeschiari et al.l 120091) on 10,000 synthetic data set realisations. The results are given in 
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Table 01 and the most recent cycles inclnding the Weihai data are plotted in Figure [H 
The rms about the £t for the Weihai velocities is 7.3 ms“^, comparable to the previously 
published da t a, an d consistent with the stellar oscillation amplitude of 5-6 ms“^ found by 


Hatzes et al.l (120121) . This result demonstrates that the Weihai data acquisition and Doppler 


velocity extraction techniques are robust. 


4. Summary and Conclusions 

Close-in planets (a < 0.5 AU, P<100 days) are rare around giant-branch stars. We 
have begun an observational program at the Im telescope located at Weihai Obervatory 
of Shandong University, using its stabilised echelle spectrograph to obtain precise Doppler 
velocity measurements of a sample of bright giant stars. We aim to observe these stars with 
as high a cadence as possible to search for these “missing” planets. Preliminary results from 
this progam give excellent results for the known planet-host Beta Gem (HD 62509), and 
conhrm our ability to obtain precise velocities with this new instrument. 


The Weihai iodine cell is provided by Bun’ei Sato and the Okayama Astrophysical 
Observatory. RW acknowledges support from UNSW Faculty Research Grants. S. M. Hu 
would like to thank the support by the National Natural Science Foundation of Ghina and 
Ghinese Academy of Sciences joint fund on astronomy under grant No. U1331102, by the 
National Natural Science Foundation of Ghina under grant No. 11333002 and by Sino- 
German Science foundation under project No. GZ788. M.E. is supported by the National 
Science Foundation through grant AST-1313075. 
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HIP 

3031 

4422 

4906 

5586 

6411 

7294 

9884 

13061 

14668 

14838 

19038 

20252 

20455 

20877 

20885 

20889 

22957 

24822 

26366 

26885 

27483 

28358 

29696 

37740 

37826 

39424 

42527 

42911 

47029 

51808 

53229 

58948 

59847 

60172 

62763 

63608 

72125 

73620 

74666 

75458 

77070 

77655 
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Table 1. Weihai target list 


HD RA Dec 


3546 

00 

38 

33.50 

+29 

18 

44.5 

5395 

00 

56 

40.01 

+59 

10 

52.2 

6186 

01 

02 

56.66 

+07 

53 

24.3 

7106 

01 

11 

39.59 

+30 

05 

23.0 

8207 

01 

22 

20.39 

+45 

31 

43.5 

9408 

01 

33 

55.93 

+59 

13 

55.5 

12929 

02 

07 

10.29 

+23 

27 

46.0 

17361 

02 

47 

54.44 

+29 

14 

50.7 

19476 

03 

09 

29.63 

+44 

51 

28.4 

19787 

03 

11 

37.67 

+ 19 

43 

36.1 

25604 

04 

04 

41.66 

+22 

04 

55.4 

27348 

04 

20 

24.66 

+34 

34 

00.3 

27697 

04 

22 

56.03 

+ 17 

32 

33.3 

28292 

04 

28 

26.37 

+ 16 

21 

34.7 

28307 

04 

28 

34.43 

+ 15 

57 

44.0 

28305 

04 

28 

36.93 

+ 19 

10 

49.9 

31421 

04 

56 

22.32 

+ 13 

30 

52.5 

34559 

05 

19 

16.59 

+22 

05 

48.1 

37160 

05 

36 

54.33 

+09 

17 

29.1 

37984 

05 

42 

28.66 

+01 

28 

28.8 

38656 

05 

49 

10.46 

+39 

10 

52.1 

40035 

05 

59 

31.55 

+54 

17 

05.9 

43039 

06 

15 

22.74 

+29 

29 

55.4 

62345 

07 

44 

26.87 

+24 

23 

53.3 

62509 

07 

45 

19.36 

+28 

01 

34.7 

66216 

08 

03 

31.10 

+27 

47 

39.9 

73108 

08 

40 

12.90 

+64 

19 

40.3 

74442 

08 

44 

41.11 

+ 18 

09 

17.5 

82741 

09 

35 

03.85 

+39 

37 

17.2 

91190 

10 

35 

05.59 

+75 

42 

46.7 

94264 

10 

53 

18.64 

+34 

12 

56.0 

104979 

12 

05 

12.67 

+08 

43 

58.2 

106714 

12 

16 

20.56 

+23 

56 

43.5 

107328 

12 

20 

21.15 

+03 

18 

45.8 

111812 

12 

51 

41.93 

+27 

32 

26.6 

113226 

13 

02 

10.76 

+10 

57 

32.8 

129972 

14 

45 

14.50 

+16 

57 

51.9 

133165 

15 

02 

54.07 

+02 

05 

28.6 

135722 

15 

15 

30.10 

+33 

18 

54.4 

137759 

15 

24 

55.78 

+58 

57 

57.7 

140573 

15 

44 

16.00 

+06 

25 

31.9 

142091 

15 

51 

13.94 

+35 

39 

29.6 


V 

Teff (K) 

log 9 (cgs) 

4.34 

5102 

2.80 

4.62 

4875 

2.7 

4.27 

4955 

2.68 

4.51 

4748 

2.62 

4.87 

4656 

2.8 

4.68 

4900 

2.67 

2.01 

4498 

2.4 

4.52 

4727 

2.67 

3.79 

5022 

3.12 

4.35 

4732 

2.7 

4.36 

4699 

2.6 

4.93 

5073 

3.10 

3.77 

5058 

3.02 

4.96 

4529 

2.5 

3.84 

4955 

2.9 

3.53 

4797 

2.6 

4.06 

4498 

2.4 

4.96 

5096 

3.17 

4.09 

4898 

2.9 

4.90 

4508 

2.2 

4.51 

5021 

2.90 

3.72 

4911 

2.86 

4.32 

4732 

2.7 

3.57 

5101 

3.04 

1.16 

4998 

3.13 

4.94 

4560 

2.5 

4.59 

4564 

2.28 

3.94 

4763 

2.67 

4.81 

4934 

2.81 

4.86 

4890 

3.07 

3.79 

4677 

2.8 

4.12 

4996 

2.86 

4.93 

4887 

2.6 

4.97 

4514 

1.94 

4.93 

5623 

2.9 

2.85 

5145 

3.12 

4.60 

4887 

2.7 

4.39 

4825 

2.79 

3.46 

4963 

2.73 

3.29 

4477 

2.5 

2.63 

4498 

2.5 

4.79 

4764 

3.0 


[Fe/H] Reference 

-0.56 (200T) 

-0.51 ^assarott^^al^ ^00^ 

-0.24 ^uc^^^leit^ (^007) 

0.01 ^uc^^^leitCT (200T) 

0.03 ^assarott^^al^ ^00^ 

-0.22 Luck fc Heiter (2007) 

-0.25 ^assarott^^al^ (^00^ 

0.07 Luck h Heiter ^2007) 

0.19 ^uc^^^leit^ (^007) 

-0.03 ^assarott^^al^ ^00^ 

0.01 ^assarott^^al^ ^00^ 

0.09 ^uc^^^leitCT (^007) 

0.19 ^ucl^^^leitCT (^007) 

-0.17 ^assarott^^al^ (^00^ 

0.04 ^assarott^^al^ ^00^ 

0.04 ^assarotti^^al^ ^00^ 

-0.26 ^assarott^^al^ ^00^ 

0.10 ^uc^^^leitCT (^007) 

-0.63 ^assarott^^al^ ^00^ 

-0.55 ^assarott^^al^ (^00^ 

-0.09 ^uc^^^leitCT (^007) 

-0.01 ^uc^^^leit^ (^007) 

-0.33 ^assarott^^al^ (^00^ 

0.03 ^uc^^^leitCT (200T) 

0.17 ^uc^^^leitCT (^007) 

0.03 ^assarott^^al^ (^00^ 

-0.16 ^ucl^^^leitCT (200T) 

0.01 ^uc^^^leit^ (^007) 

-0.12 ^uc^^^leitCT (200T) 

-0.15 ^cWilliam (199^ 
-0.20 ^assarott^^al^ (^00^ 

-0.33 ^uc^^^leit^ (^007) 

-0.24 ^assarott^^al^ (^00^ 

-0.46 ^uc^^^leit^ (^007) 

0.01 ^assarott^^al^ ^00^ 

0.13 ^uc^^^leitCT (200T) 

-0.10 ^assarott^^al^ ^00^ 

-0.13 ^uc^^^leit^ Q007) 

-0.30 ^ucl^^^leitCT (200T) 

0.03 ^assarott^^al^ ^00^ 

0.03 ^assarott^^al^ ^00^ 

-0.04 ^assarott^^al^ (^00^ 
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Table 2. Known Planet Hosts in the WES Sample 


Planet 

Period (days) 

M sin i (Mjup) 

a (AU) 

Reference 

a Ari (HD 12929) 

380.0±0.3 

1.72±0.19 

1.130±0.062 

Lee et al. f2011^ 

e Tau (HD 28305) 

595±5 

7.7±0.3 

1.936±0.034 

Sato et al. ('20071 

/3 Gem (HD 62509) 

589.6±0.8 

2.76±0.14 

1.757±0.029 

Hatzes Cochran f 19931 

4 UMa (HD 73108) 

269±2 

7.1±0.6 

0.877±0.036 

Dollineer et al. 120071 

i Dra (HD 137759) 

511.098±0.089 

9.3±0.9 

1.31±0.06 

Frink et al. ('20021 

K CrB (HD 142091) 

1300±15 

1.97±0.12 

2.72±0.05 

Johnson et al. ('20081 


Table 3. Weihai radial velocities for HD 62509 


BJD-2400000 

Velocity (ms 

Uncertainty (ms 

56293.24524 

-9.0 

7.7 

56293.25392 

-6.1 

7.9 

56293.26226 

-0.3 

7.9 

56293.27088 

-7.5 

7.4 

56316.09065 

-11.9 

8.4 

56316.10066 

-2.0 

7.9 

56316.11256 

-13.7 

8.2 

56316.12067 

-2.4 

8.6 

56316.12989 

-12.9 

8.2 

56317.12107 

-6.8 

7.2 

56317.13220 

-3.0 

7.8 

56320.13640 

-16.6 

8.0 

56320.14513 

-12.8 

8.4 

56320.15348 

-11.7 

8.1 

56320.16399 

-11.6 

8.5 

56320.17287 

2.6 

8.7 

56617.27237 

73.8 

9.2 

56617.28150 

72.2 

8.9 

56617.28955 

72.6 

9.0 

56617.29755 

74.0 

9.1 

56617.30750 

83.4 

9.1 

56617.31558 

84.1 

9.9 

56617.32360 

81.1 

8.7 

56617.33160 

82.5 

8.9 

56617.34040 

79.5 

9.4 

56617.34844 

79.0 

9.4 

56724.03493 

43.2 

10.2 

56749.01689 

46.8 

8.1 

56749.04009 

48.4 

8.7 

56749.07138 

43.2 

8.1 

56763.99622 

40.0 

7.8 

56764.02025 

46.0 

8.3 

56781.98304 

27.1 

8.6 

56781.99446 

31.3 

8.6 

56782.03382 

2.5 

8.4 

56786.99866 

32.3 

8.2 

56792.98783 

13.2 

9.5 

56792.99932 

37.6 

10.8 
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JD-2400000 


Fig. 1.— WES data for Beta Gem: 38 observations from 2013 Dec 31 to 2014 May 15. Our 
results are in excellent agreement with the published data, now spanning 33.5 years and 
conhrming the consistency of the planet’s orbital parameters. 















Table 4. Updated Parameters for HD 62509b 


Parameter 

Value 

Period 

591.2±0.76 days 

Eccentricity 

0.071±0.028 

UJ 

263±19 degrees 

K (ms“^) 

44.1±1.2 ms-l 

To 

2444036.3±32.5 JD 

m sin i 

2.80±0.08 Mjup 

a 

1.7113±0.0015 AU 

RMS - McD cs21 

14.9 ms“^ 

RMS - McD phase 3 

13.7ms-l 

RMS - TOPS 

12.0ms-l 

RMS - CFHT 

19.2ms-l 

RMS - DAO 

37.2 ms-l 

RMS - McD 2.1m 

23.2 ms-l 

RMS - Lick 

9.2ms-l 

RMS - Weihai 

7.3ms-l 







